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Abstract. We present the CCD Cousins R band photometric observations of the 
afterglow of GRB 041006. The multiband afterglow evolution is modelled using an 
underlying 'hard' electron energy spectrum with a pi ~ 1.3. The burst appears to be of 
very low energy {E ~ lO'** ergs) confined to a narrow cone of opening angle 9 ~ 2.3°. 
The associated supernova is compared with SN1998bw and is found to be brighter. 
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1. Introduction 

GRB 041006 was detected by the High Energy Transient Explorer (HETE II) Fregate and WXM 
instruments at 12:18:08 UT on 06 October 2004 (Galassi et al. 2004). The refined error circle 
provided by the HETE II ground analysis was of radius 5 arc minutes centered at a = 00^54™53'', 
5 = +01° 12' 04" (J2000). The fluence of GRB 041006 in 2-30 keV and 30-400 keV were 5 x 10"^ 
erg/cm^ and 7 x 10~^ erg/cm^ respectively which classifies this as an "x-ray rich GRB" . The light 
curve shape and the spectral characteristics of GRB 041006 were similar to GRB 030329 though it 
was 20 times fainter than GRB 030329. A fading x-ray counterpart was reported by Butler et al. 
(2004) and the optical afterglow located at a = 00h54'"50M7, S = +01°14'o7" (J2000) was found 
by Da Costa, Noel and Price (2004). From low resolution spectroscopy performed on October 
7.10 UT Fugazza et al. (2004) derived a redshift of z = 0.712 for this burst. This value was 
later confirmed by Price et al. (2004). The optical afterglow showed a 'supernova bump' around 
5.0 days. Adopting the fluence reported by Galassi et al. (2004) and the estimated redshift, 
Fugazza et al. (2004) derive the isotropic equivalent gamma— ray energy of the burst to be Eiso = 
1.5 X 10^^ erg. 

In this paper we present Cousins R filter optical observations of the afterglow of GRB 041006 
carried out during the first night of the GRB trigger and the modelling of the data along with the 
multiband observations available so far in the literature. For a reliable determination of the optical 
afterglow magnitudes we imaged the field of GRB 041006 along with RU 152 standard region of 
Landolt (1992) and calibrated a total of 10 stars in the field. The optical observations and the data 
reduction is discussed in section 2 and the modelling assuming a 'double slope' electron energy 
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spectrum is discussed in section 3. Section 4 explains briefly about the associated supernova and 
its comparison with SN1998bw. A brief discussion is presented in section 5 and the conclusions 
follow in section 6. 



2. Observations and Data Reduction 

The observations in Cousins R filter of the optical afterglow of GRB 041006 were carried out on 
06 October 2004 using the 104-cm Sampurnanand Telescope of ARIES, Nainital. One pixel of the 
2048 X 2048 pixeP size chip corresponds to a square of 0".38 side, and the entire chip covers a 
field of 13' X 13' on the sky. The gain and read out noise of the CCD camera are 10 e~/ADU and 
5.3 e~ respectively. 

We observed the Landolt (1992) standard RU 152 and the OA field in VRI filters on 08 October 
2004 for photometric calibration during good photometric sky conditions. Several bias and twilight 
flat frames were observed with the CCD camera to calibrate the images using standard techniques. 

ESO MIDAS, NOAO IRAF and DAOPHOT softwares were used to process the images. The 
values of the atmospheric extinction on the night of 08/09 October 2004 determined from the 
observations of RU 152 bright stars arc 0.19, 0.15 and 0.11 mag in V, R and I filters respectively. 
The 7 standard stars in the RU 152 region cover a range of -0.06 < (T^ — ii) < 0.47 in color and 
of 11.08 < y < 15.02 in brightness. This gives the following transformation equations 

v^^a = V - (0.11 ± 0.01)(y -R) + (4.29 ± 0.004) 
r^cd = R- (0.12 ± 0.01)(U - R) + (4.20 ± 0.004) 
i^^^ = I- (0.14 ± 0.02)(i? - J) + (4.69 ± 0.007) 

where V, R, I are standard magnitudes and Vccd^ Tccd^ iced represent the instrumental magnitudes 
normalized for 1 s of exposure time and corrected for atmospheric extinction. The color coefficients, 
zero-points and errors in them were determined by fitting least square linear regressions to the 
data points. Using these transformations, VRI photometric magnitudes of 10 secondary stars were 
determined in GRB 041006 field and their values are listed in Table 1. The {X,Y) CCD pixel 
coordinates are converted to a2ooo, <^20oo values using the astrometric positions given by Hcnden 
(2004). The 10 secondary stars in the GRB 041006 field were observed 2, 4 and 2 times in V,R 
and / filters respectively. These stars have internal photometric accuracy better than 0.01 mag. 
The zero-point differences on comparison between our photometry and that of Hendcn's (2004) are 
0.058 ± 0.043, 0.062 ± 0.045 and 0.026 ± 0.041 in V, R and / filters respectively. These differences 
are based on the comparison of the 10 secondary stars in the GRB 041006 field. 

Our first observation was made at ~ 0.24 days after the burst. Several short exposures, with 
exposure time varying from 5 min to 25 min, were taken to image the OA of GRB 041006. To 
improve the signal-to- noise ratio of the OA, the data have been binned in 2x2 pixel^. Profile- 
fitting magnitudes were determined from these images using DAOPHOT II software. The profile 
magnitudes were converted to aperture (about 4 arcsec) magnitudes using aperture growth curve 
determined from well isolated secondary stars. They were differentially calibrated with respect to 
the secondary stars listed in Table 1. The magnitudes derived in this way are given in Table 2. 

The published photometric measurements by Kinugasa ct al. (2004a), Ferrcro ct al. (2004), 
Monfardini et al. (2004), Fynbo et al. (2004), Fugazza et al. (2004), Greco et al. (2004), D'Avanzo 
et al. (2004), Garnavich ct al. (2004), Balman et al. (2004), Covino et al. (2004), Bikmaev et 
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Table 1. The identification number (ID), (a, 5) for epoch 2000, standard V, {V - R) and {R - I) 
photometric magnitudes of the stars in the GRB 041006 region are given. Number of observations taken 
in V, R and I filters are 2, 4 and 2 respectively. 



ID 


ci;20oo 


6 


2000 


V 


V-R 


R-I 




(h m 


s) 


(de. 


g m s) 


(mag) 


(mag) 


(mag) 


1 


00 54 


47 


01 


14 05 


17.99 


0.19 


0.27 


2 


00 54 


37 


01 


14 02 


13.67 


0.33 


0.39 


3 


00 55 


01 


01 


11 25 


13.98 


0.43 


0.40 


4 


00 54 


38 


01 


12 59 


17.36 


0.57 


0.54 


5 


00 54 


54 


01 


15 47 


16.96 


0.59 


0.56 


6 


00 54 


53 


01 


17 47 


15.19 


0.58 


0.53 


7 


00 54 


55 


01 


16 58 


15.21 


0.42 


0.40 


8 


00 54 


57 


01 


19 31 


14.38 


0.60 


0.53 


9 


00 55 


05 


01 


17 07 


16.21 


0.38 


0.36 


10 


00 55 


01 


01 


15 39 


15.47 


0.47 


0.47 



Table 2. CCD Cousins R band optical observations of the GRB 041006 afterglow using the 104-cm 
Sampurnanand Telescope at ARIES, Nainital. 



Date (UT) 


Magnitude 


Exposure Time 


Passband 


2004 October 


(mag) 


(s) 




6.7572 


19.80 ± 0.056 


300 


R 


6.7666 


19.74 ± 0.034 


600 


R 


6.7763 


19.91 ± 0.043 


600 


R 


6.7867 


19.93 ± 0.053 


800 


R 


6.8151 


20.02 ± 0.047 


1200 


R 


6.8305 


20.06 ± 0.059 


1200 


R 


6.8472 


20.04 ± 0.047 


1200 


R 


6.8624 


20.06 ± 0.059 


1200 


R 


6.8772 


20.23 ± 0.062 


1200 


R 


6.9101 


20.22 ± 0.053 


1500 


R 



3. Multiband Modelling 



The afterglow flux fi,{t) observed at a given frequency v at time t after the burst is usually described 
by a powerlaw (/i/(t) oc u~^t~"), with a and /3 assuming different values in different spectral and 
dynamical regimes of evolution (See Piran 2004 for a recent review). The optical afterglow of 
GRB041006 exhibited a relatively flat decay ax ~ 0.5 in the early phase which steepened to 
a2 ~ 1.3 after a break around 0.14 day (Stanek et al. 2004, Soderberg et al. 2005). If this break 
is due to the lateral spreading of a conical fireball, the indices indicate a relatively hard (p < 2) 
energy spectrum for the synchrotron emitting electrons (Bhattacharya 2001). At 0.058 day the 
intrinsic spectral index of the afterglow was 0.685 ±0.13 between B and V bands (Da Costa & Noel 
2004). At 1.74 day, the index became 1.0 ± 0.18 between B and R bands (Garnavich ct al. 2004). 
But this is only a marginal steepening, seen within one sigma error. The x-ray spectral index of 
0.9 ± 0.2 measured at 1.23 days (Butler et al. 2004) is consistent with the /3 observed within the 
optical bands. Hence, the reported spectral indices in both optical and x-ray bands support the 
inference of a hard electron energy spectrum. Butler et al. (2004) report that the x-ray lightcurve 
follows a single power law of slope — 1.0 ± 0.1 during 0.7 to 1.7 days. This indicates that the x-ray 
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Figure 1. A 12' x 12' R band image of GRB 041006 optical afterglow. The position of the optical 
afterglow is marked with a circle. The 10 secondary stars used for calibration are also marked. North is 
up and East is to the left. 

from the optical and x-ray observations, upper limits to the afterglow flux are available in near IR 
(Kinugasa et al. 2004b), millimeter (Barnard et al. 2004a, 2004b) and also in radio (Soderberg et 
al. 2004) wavelengths. 

We attempted modelling the multi wavelength afterglow evolution with an underlying 'double 
slope' electron energy spectrum. The electrons arc assumed to be distributed in a flatter spectrum 
(pi < 2) until a cutoff energy jiirieC^, and thereafter towards the higher energy end, the distribution 
steepens to p2 > 2 (see Panaitescu and Kumar (2001) for a discussion on such double slope electron 
energy distribution), yi is termed as the 'injection break'. The evolution of ji significantly affects 
the afterglow light curve. Following Bhattacharya (2001), we parametrize the evolution of the 
injection break as 7i = ^F' where F is the bulk Lorentz factor of the shock and ^ is the constant 
of proportionality. 

We assume the cooling frequency (fc) to be below the optical bands to satisfy a of 0.5 and (3 
in the range of 0.6 — 0.7 simultaneously. There is no signature of a steepening seen at the higher 
energy end of the spectrum from the available observations. The upper cut off of the electron 
energy spectrum may lie at energies higher than what the observations could probe. Hence we 
place the 'injection break frequency', i>i, corresponding to the electron Lorentz factor 7^, above 
the x-ray band. We compute the spectral evolution of the afterglow with these basic assumptions. 
A broken power-law model for the underlying spectrum and evolution is used (Sari et al. 1998 ), 
with Band- type smoothcning (Band et al. 1993) applied to spectral and temporal breaks. After 
excluding the NIR/mm/radio upper- limits and optical data beyond 5 days (where a supernova 
bump prominently appears) from the fit, the model converged to a minimum Xdof 3.4. An 
additional extinction other than that due to Galactic dust {Ay of 0.07 mag, estimated for this 
direction using the data of Schlegel et al. 1989) is required by the model, which could result from 
the dust distribution in the host galaxy of the GRB. We obtained Ay in the range of 0.03— 0.16 mag 
in the host rest frame, assuming a Galaxy type extinction law for the dust distribution (Cardelli 
et al. 1989). According to Butler et al. (2004) the x-ray spectrum is well fit by an absorbed 
power-law with the neutral hydrogen column density (Nh) greater than that expected from the 
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Table 3. Model Parameters 



Electron Distribution 
Parameters 


Spectral 
Parameters 
(At 0.1 days) 


Pi 


1.29-1.32 




(1.2 - 3.0) X 10^^ Hz 




> 2.2 


J^c 


(1.0 - 2.0) X 10^4 Hz 


q 


1.0 


Vi 


< 2.4 X 10^" Hz 




0.17 - 0.24 day 


fm 


(0.37 - 0.49) mJy 


Others 


E{B-v) (host) 0.01 - 0.05 mag 







Table 4. Physical Parameters 



n (assumed) 


1.0 atom/cc 




~ 1.0 




0.07-0.14 




> 2.0 X lO'^ 


E 


(2.0 - 4.0) X 10^^ ergs 




1.7° - 2.8° 


Etot 


(1.4-3.4) X 10*** ergs 



than the empirical value obtained for the ISM of our galaxy (Predehl & Schmitt, 1994). This 
might imply a different type of dust in the GRB host, and is not a surprising fact if GRBs are to 
occur at star forming regions, where dust depletion is highly probable. 

In determining the physical parameters, the total number of unknowns were five (isotropic 
equivalent energy release -Eiso, number density of the ambient medium n, fractional energy content 
in electrons Cg, that in the magnetic field cb and ^, the constant of proportionality for 7^. However 
the total number of parameters inferred from the observations are four since the self absorption 
frequency Va is not well determined. Wo hence assiimed a range of values for n, and obtained the 
remaining four physical parameters for the afterglow. Since only lower limit is available for Vi and 
P2, C is not well constrained from the upper end. This will also affect the precise determination 
of eg, but the best fit provides a value close to unity for e^. Table 3 & 4 lists the fit parameters 
and the derived physical parameters for a nominal value of n = 1.0 respectively, from the best fit 
model. A decrease in the value of n to 0.05 atom/cc increases Eiso to 8 x 10^^ ergs, and cb to 0.2. 
By increasing the number density to 50 atom/cc, the isotropic energy decreased to 10''''^ ergs and 
cb decreased to 3. x 10~^. We obtain a range of (1.0 — 7.0) x 10^^ ergs for the total energy release 
corresponding to a range of 0.05 — 50.0 atom/cc for n. We used z = 0.716 (Price et al. 2004) and 
a corresponding luminosity distance of 4.72 Gpc (assuming a cosmology of Hq = 65 km/sec/Mpc, 

= 0.3 and f^A = 0.7). We have displayed the observations along with the model in figure 2. 

Since available observations do not cover the 0.12 day break in various bands, its achromatic 
nature as expected for a jet break, cannot be confirmed. As a mild evolution in the optical spectral 
index between 0.06 days to 1.7 days has also been reported, we examined the possibility of the 
temporal steepening arising from the passage of a spectral break through R-band. The observed 
change in a is ~ 0.8 and the maximum change that could be attributed to /? is 0.3. Passage of 
the cooling frequency would predict Aa = 0.25 and A/3 = 0.5, and hence can not reproduce the 
observed a and /3. The only other possibility is the passage of the injection break through R-band 
around 0.1 day, with pi being ~ 1.3 and p2 ~ 2.3. A model with this assumption can reproduce 
the optical observations and the NIR/mm/radio upper limits very well but it underpredicts the 
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Figure 2. GRB 041006 afterglow observations in different bands along with the model. Filled circles 
represent our data and open circles represent the data from literature. Light curve in different bands are 
shifted as indicated. Inset : X-ray observations and radio upper limits along with the model. 

p < 2, we calculated the inverse compton emission for this model, and found that it is negligible 
at the x-ray frequencies. 

Wc find from our analysis that the best possible explanation for the multi wavelength obser- 
vations for this afterglow is the synchrotron emission from a coUimated outflow with electrons 
distributed in a hard energy spectrum. 



4. The associated supernova 

The afterglow optical flux decay has shown deviation from a power law behaviour after ~ 0.5 
days and exhibited a bump around 25 days. Stanek et al. (2004) attributes this behaviour to 
the emission from an associated supernova. We subtracted the modelled power-law flux and an 
assumed host galaxy flux from the observed emission to estimate the supernova contribution. 
There is no measurement available for the host galaxy flux and we are left with only assumptions. 
Stanek et al. (2004) find that the supernova associated with GRB 041006 peaks at a later (1.35 
times) time after explosion and is also brighter in comparison to SN1998bw redshifted to z = 
0.716. Soderberg et al. (2005) from the late HST observations (mainly I band) conclude that the 
supernova is ~ 0.3 magnitudes brighter than a k-corrected SN1998bw whilst no additional stretch 
in time is required. They also point out the possibility of late time flux getting contaminated by 
the host and also by a nearby galaxy. 

The major difference our afterglow model has, in this context, is the presence of the host 
extinction (rest frame Ay = 0.1 mag in the best fit model). We used the late afterglow magnitudes 
of Stanek et al. (2004) in R-band and Soderberg et al. (2005) in I band to obtain the residual flux, 
which could be attributed to the associated supernova after subtracting the host flux. We assumed 
four free parameters (host galaxy flux in both the bands. Am, the shift in magnitude required 
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Figure 3. Comparison between tiie associated supernova in R (squares) and I (triangles) bands witli tlie 
redshifted k-corrected SN1998bw (line) after applying a Am and 5 correction (see text). The parameters 
used in obtaining tho data sets are Am = —0.65, 5 = 1.05, R(host) = 26.2 and I(host) = 28.9. The curves 

have been shifted as indicated. 

to —0.75, S between 1.05 and 1.1, R(host) in the range of 26.2 to 28.8 mag and I(host) between 
26.8 and 28.9 mag (figure 3). Prom the figure, it appears that the I band data fits better with the 
template compared to the R band. Wc notice that a different value of Am and S would be required 
to fit the R band better, which suggests that the associated SN may have a different spectrum 
compared to that of SN1998bw. This could also be the reason for the different conclusions arrived 
at by Stanek et al. (2004) and Soderberg et al. (2005) about the associated SN. 



5. Discussion 



The decay index of ~ 1.3 seen in GRB 041006 afterglow lightcurve, post the one day break, 
suggests that the underlying electron energy spectrum could be hard (p ^ 1.3). A different model, 
proposed by Granot, Ramirez- Ruiz & Perna (2005), attributes this behaviour to the jet being 
viewed off-axis. In this model, the ~ 0.1 day break is caused by the lorentz factor falling below 
the inverse of the viewing angle, there is no lateral expansion of the jet and the derived electron 
energy distribution is steeper (p ~ 2.2). They invoke the presence of a stellar wind driven ambient 
density profile to explain the reported x-ray decay index {ux 1.1) which is flatter in comparison 
to the simultaneous optical decay index {uo ^ 1.3). Granot et al. (2005) derive a kinetic energy 
of 10^^ ergs for this explosion. 

Since we do not find the x-ray decay index to be flat (sec section 3.), wc used a normal ISM 
density profile for the ambient medium. The kinetic energy derived for the explosion is rather 
small (10** ergs). With the available set of data, it is impossible to discriminate between the 
two models. However, both off-axis jets and wind driven density profiles, although realistic, are 
not inferred commonly from the afterglow models. For this afterglow, in fact, there is no need 
of going beyond the simplistic assumption of an on-axis jet ploughing through a constant density 
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Butler et al. (2004) extrapolate the observed optical spectrum, which has an index /3o of 
— 1.0 ± 0.1 to the x-ray frequencies and find that the observed x-ray flux is underpredicted. They 
suggest that a possible explanation for this and also for the smaller could be the presence of 
considerable amount of inverse compton emission in the x-ray bands. They also find that the 
spectral index within the x-ray band (/3x ~ 0.7) is smaller than that within the optical band 
(/3o « —1.1 ±0.1). Since in our model, wc have a hard electron energy spectrum, (3 is close to 
0.7, and the presence of extinction in optical bands due to the dust column of the host galaxy, 
explains the observed higher optical (3. Hence our model does not underpredict the observe x-ray 
flux, eventhough, the inverse compton emission predicted by the model is negligible at the x-ray 
frequencies. 

6. Conclusions 

• In this paper, we present the Rc band optical observations of the afterglow. 

• We also show that multi band observations can be well reproduced by synchrotron emission 
from a hard electron energy spectrum and an early jet break {tj = 0.2 day). 

• The model requires additional extinction due to the dust distribution in the host. The E(B- 
V) derived by us is two orders of magnitude lesser than what one would predict from the Nh 

column density inferred by Butler et al. (2004) assuming a Galaxy type gas to dust ratio. 
This difference might indicate a considerable difference in the characteristics of dust column 
of the host. 

• The kinetic energy of the burst inferred from the model is low (^ 10^* ergs), and is an order 
of magnitude smaller than the energy emitted in gamma rays. 

• Comparison of the afterglow subtracted optical flux beyond 5.0 days with a redshifted 

k-corrccted SN1998bw light curve shows that the associated supernova is brighter than 
SN1998bw, and has a longer rise time and might have a different spectrum than SN1998bw. 
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